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1. 


Introduction: 


Many  important  problems  of  engineering  and  management  are  of  a  form 
which  could  be  represented  as  geometric  programs  except  that  the 
functional  to  be  minimized  as  well  as  the  constraints  are  not  confined 

If 

to  "posynomials"  in  that  some  of  the  coefficients  are  negative.  The 

resulting  problem  thus  may  not,  in  general,  be  transformed  to  an 

2/ 

equivalent  convex  programming  problem.  To  date  the  only  general 
method  for  obtaining  global  optima  to  (necessarily  non-convex)  problems 

3/ 

with  multiple  local  optima  is  Gomory's  integer  programming  method. 

We  are  herewith  proposing  an  approximate  method  for  another  class  cf 
problems  with  multiple  local  optima--viz.  ,  extensions  of  geometric 
programming  in  which  some  of  the  coefficients  are  negative.  This  method 
provides,  at  each  stage,  a  convex  approximant  which,  a  fortiori,  provides 
the  duality  relations  that  are  needed  for  many  purposes.  This  is  in 

if 

contrast  to  other  approaches  which  either  lose  these  duality  relations 

5/  ‘ 

or  else  restrict  the  applications  to  special  situations.  More  specifically, 


U  Cf.  [  7  ]  for  definitions  of  this  and  other  terminology  in  geometric 
programming. 

2/  Cf.  ,  e.  g.  ,  the  exponential  transformations  used  in  [  3  ]  and  [  4  ]. 

3/  See  [  8  ]  and  [  9  ]  for  Gomory's  original  articles.  See  also  [  2] 
and  [  6]  for  further  discussion  and  development. 

4/  Cf.  ,  e.g.  ,  [  10]. 


5/  The  constraints  in  [3]  and  [5]  ,  for  instance,  were  arranged  so  that 
they  could  always  be  treated  in  a  manner  which  did  not  preclude 
access  to  the  indicated  duality.  Other  possibilities  are  also 
present,  however,  as  witness  some  of  the  examples,  treated  in  [7]. 
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the  method  that  we  shall  describe  here  Is  conceived  in  the  same  spirit 
as  previous  suggestions  we  have  made  as  a  result  of  other  research  we  have 


1/ 

conducted  to  extend  the  boundaries  of  ordirary  linear  programming. 


2 .  Formulation  and  Development  of  tihe  Convex  Approximant: 

2/ 

Consider  the  following  problem 


(i.3) 


mln  80  -  g0 


subject  to 


“  gt  <  1 ,  i  =  l ,  .  •  •  ,  m 


where  the  g^i  g^  are  posynomials  in 


(1.2) 


t  *  (t^ ,  ....  t^)  . 


X  •  6  •  t 
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je^  J  keK1 
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1 J  Cf.  ,  e.g.  ,  [  1  ]  and  [  5  ]. 

2/  To  abbreviate  this  part  of  the  development,  it  is  assumed  that  all 
conditions  for  existence  and  attainment  of  the  indicated  minima  are 
fulfilled,  Cf.  [  7  ]  for  a  rigorous  treatment  of  the  relevant 
necessary  and  sufficient  conditions  in  complete  detail. 


I 
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Note  that  the  above  problem  i8  a  generalization  of  ordinary  geometric 
programming  in  that  the  constraints  and  the  functional  are  not  confined 
to  posynomials. 

3.  Formulation  of  Approximants : 

Each  one-term  posynomial  P^(t)  preceding  expressions  mav 

be  replaced  by  a  single  variable  y  subject  to 


(2.1) 


ylj  -  Plj^ 


or 


(2.2) 


yij  Cpij ( 1  <  1 


which  is  the  same  as 


(2.3) 


'ij 


-b1J  -bij 

,t  1  .  .  .  t  n 

LI  n 


<  1  . 


The  resulting  problem  in  t  and  the  y^  is  equivalent  to  (1.1). 

Next,  let  us  suppose  that  the  range  of  each  y^j  relevant  to  tha 
optimization  may  be  represented  by 


(3) 


o  <  <  ytJ  <  «tj 


We  then  introduce  k  >  U  and  consider  the  function 


(4) 


W  :  ku  -yij 


as  diagrammed  below.  Evidently  over  the  interval  (L^  ,  U^j)  the  linear 
function  (4)  is  positive  and  bounded  above  and  below.  It  may  thus  be 


i 
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approximated  by  a  polynomial 


(5) 


W 


To  the  degree  of  approximation  thus  rendered--e.  g.  ,  approximation 
of  the  linear  function  by  posynomials — the  original  problem  (1.1)  is  now 
replaced  by 
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(6) 


m 


min  g+(t>+^  qoJ  <yoJ) 


subject  to 


m 


gi(t)  +  ^  qu  (yij>  ^ 1  +  j  ku 


[p±j  (t)l  1  ytJ  <i 

yU  "ij  *  1 

yij  Lij  -  1 


t  >  0 


This  problem  may  evidently  be  transformed  (e.g.,  by  the  exponential 

1/ 

transformation)  into  a  convex  programming  problem.  We  therefore 
call  It  a  convex  approximant  of  the  original  problem.  It  therefore 
follows  that  it  has  only  one  local  (=  global)  optimum  value. 

Note  in  particular  that  each  convex  approximant  has  an  associated 
dual  problem.  Thus  a  dual  evaluator  is  available  for  each  constraint. 
Those  that  refer  to  the  constraints  indicate  possible  directions 

of  improvement  if  these  upper  or  lower  bounds  are  tight.  The  dual 
evaluator  is,  of  course,  equal  to  zero  when  These  bounds  are  slack.  The 
approximation  can  thus  be  improved  in  the  neighborhood  of  any  already 
attained  optimum  by,  e.g.  ,  reducing  the  range  of  the  slack  and  , 

thereby  enabling  one  to  make  an  improved  posynomial  fit  in  the  next 


1/  See  [  3  ]  and  [4  ]. 


-  — y 
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convex  programming  approximant.  Similarly,  the  interval  may  be 
reduced  and  translated  in  the  direction  indicated  by  the  non-zero  dual 
evaluator  for  the  tight  ,  L^j  constraints. 

Thus,  sequentially,  the  convex  approximant  can  be  refined.  One 
would  expect  the  global  optimum  to  be  obtained  by  this  method  in 
situations  where  the  original  problem  has  multiple  local  optima.  For, 
if  the  global  optimum  value  were  significantly  different  from  that  of 
other  local  optima,  one  would  anticipate  that  the  small  modifications 
of  the  smooth  continuous  functionc  to  equally  smooth  continuous 
approximants  would  not  significantly  alter  the  global  optimum.  Since 
the  convex  approximant  has  onl’"  one  local  (=  global)  optimum,  its  value 
should  therefore  be  close  to  the  global  optimum  value  of  the  original 
problem.  On  the  other  hand,  when  the  global  optimum  value  of  the 
original  does  not  differ  significantly  fromother  local  optimum  values, 
the  precise  optimum  obtained  matters  little  so  far  as  value  is  concerned. 
In  either  situation  therefore  one  would  expect  a  sequence  of  convex 
approximants  to  yield  a  worthwhile  result. 

3.  Conclusion: 

In  the  paper  [4  ] ,  we  showed  how  geometric  programming  could 
be  applied  to  the  determination  of  multiple  simultaneous  EOQ  (economic 
order  quantity)  formulas  under  constraints  as  well  as  to  aspects  of  the 
economic  theory  of  production  (e.g. ,  with  Cobb-Douglas  and  generalized 
SMA.C  production  functions).  Still  further  extensions  in  this  direction 
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(e.g.  ,  to  problems  of  capital  budgeting)  critically  depend  on  the 

possibility  of  dealing  with  the  presence  of  negative  coef f icients-~as 

in  (1.1)—  and  the  same  is  true  even  of  the  originally  motivated  applications 

to  engineering  designs  when,  for  instance,  scrap  values  require 

consideration.  Even  more  important,  however,  is  the  need  for  increased 

flexibility  as  when,  for  instance,  there  is  a  need  to  deal  with  problems 

where  the  natural  original  orientation  is  toward  maximization  (rather 

than  minimization)  and  where  a  restriction  to  posynomials  only  makes  it 

impossible  to  proceed  through  the  negative  of  an  associated  minimization 

1/  - -  ““  1  1  " 

problem.  A  recourse  to  the  convex  approximant  method  wouLd  then  seem  to 

be  in  order--at  least  in  these  cases  and  possibly  others  as  well. 


1/  E.g.,  as  in  ordinary  linear  programming.  Cf.  ,  e.  g.  ,  [  2  ]  or  [  6  ]. 
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